Length-dependent resistance model for a single-wall Carbon nanotube 
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The non-hnear length-dependent resistance, 7?.(Z) observed in single-wall Carbon nanotubes 
(SNTs) is explained through the recently proposed ionization energy (-E/) based Fermi-Dirac statis- 
tics (iFDS). The length here corresponds to the Carbon atoms number (A/") along the SNT. It is 
also shown that TZyily) < TZ^ifx) is associated with < Ef , which can be attributed to different 
conducting properties in their respective y and x directions, or due to chirality. 
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I. INTRODUCTION 



Enormous amount of research have been poured since 
the discovery of Carbon (C) nanotubes (CNTs) by 
lijima [l| in 1991 and consequently, CNTs have been 
successfully exploited to produce cathode ray tubes 
and nano-electronic devices 13 . Understandably, CNTs 
are believed to pave the pioneering pace for the nan- 
otechnology boom. Basically, C can be categorized into 
graphite, diamond and FuUerenes based on their bond- 
ing nature that gives rise to different electronic and struc- 
tural properties. Unexpectedly, C in all these three struc- 
tures with slight manipulations have exposed supercon- 
ductivity CNTs' electronic properties are equiv- 
alent to rolled-graphite 0, H, M, EOl which also reveal 
superconductivity in the absence of doping [ll|, [H, [III 
and concentration-dependent non-linear optical proper- 
ties. The real part of third-order non-linear susceptibil- 
ity, Re x*-^-* was found to be in the order of 10~^^ esu 
for multi-wall CNTs by Elim's group [ij]- This value is 
roughly 100 x larger than that of SNTs, which is due to 
SNT's lower C-atom concentration. The superconduct- 
ing properties of Boron-doped diamonds [l5| based on 
rcsonating-valence-bond mechanism was put forward by 
Baskaran [1, [H] whereas, the superconducting FuUerenes 
and its non-linear optical properties have been discussed 
by Cohen et al. [3 and Elim et al. [l3| respectively. 

Here, the ionization energy based Fermi-Dirac statis- 
tics (iFDS) is employed to derive the length-dependent 
resistance model, TZ{1). The derivation of iFDS and its 
applications in a wide variety of str ongl y correlated elec- 
tronic matter is given in the Refs. [3 [l^, [13, Hlj . This 
model is shown to be viable in addressing the recent 
TZ{1) observation reported by de Pablo et al. [13], An- 
driotis et al. [H] and Purewal et al. ^ in CNTs. The 
length-dependent resistance is an intrinsic property ba- 
sically because the contact resistance is independent of 
CNT's length As a consequence, the only question- 
able result is the magnitude of the resistance, not its 



length-dependent trend. However, other measurements 
namely, the temperaturc(r)-dcpendcnt electrical or heat 
conductance are strongly influenced by the contact resis- 
tance due to its own T-dependence and its large magni- 
tude, usually in the order of the CNTs resistance, which 
in turn waver the intrinsic experimental TZ{T) results. 
It is interesting to note that the resistance of a SNT is 
non-linearly proportional to the tube's length in both 
metallic and semiconducting SNTs [H, [H, [131 at any 
given T. However, the calculations carried out by Zhang 
et al. [1^ and Uryu et al. [H, ^ for metallic CNTs 
indicate that the resistance is inversely proportional to 
the length as a result of resonant tuncUing at interface. 
In this work, we do not consider hetcrostructurcs with 
resonant tunneling, but rather, on intrinsic metallic and 
semiconducting SNTs. The resistance model derived here 
are also suitable in other strongly correlated nanotubes 
that allow direct-current resistance and/or polarization 
measurements, or if the C atoms in CNTs are doped sub- 
stitutionally with different atoms. 



II. THE LENGTH-DEPENDENT RESISTANCE 
MODEL 

We start with the many-body Hamiltonian [13, [lH , 



of which. 



— VV=(i? + nr)V, 



(1) 



(2) 
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From Eq. ([2|), one can notice that the influence of 
the potential energy on the total energy has been conve- 
niently parameterized as ^. This energy function, ^ can 
be characterized in such a way that Eq is the total energy, 
E a.tT = 0. Add to that, from Eq. it is obvious that 
the magnitude of ^ is given by ^ = i?kin — ^^o + ^(r) , i?kin 
denotes the kinetic energy. Physically, ^ implies the en- 
ergy needed to overcome the potential energy that exists 
in a particular system. That is, ^ is the energy needed 



2 



to excite a particular electron to a finite distance, r, not 
necessarily r — > cx). Literally, this is exactly what we 
need to know in any condensed matter that actually or 
reasonably defines the properties of the fermions. H is 
the usual Hamilton operator, ip denotes the many-body 
eigenstate and Eq is the total energy at T = 0. The + 
sign of ±^ is for the electron (0 — ^ +00) while the — sign 
is for the hole (— cx) 0). In addition, we define the 
ionization energy in a many-atom system, ^ = EY^"^ is 
approximately proportional to Ej of an isolated atom or 
ion. We can prove the validity of Eq. ^ by means of 
constructive (existence) and/or direct proofs as given in 
Ref. [3l[. However, for an isolated atom, ^ is given by 



The corresponding total energy is 

Eo±£, = -Ekin + Vpot 

= Eq ± Ej. 



(3) 



(4) 



On the other hand, for an atom in a many-atom sys- 
tem, we can rewrite Eq. (|3]) as 



E^u 
Ei 
±E 



Eo+Vr 



pot 



many— body 



many— body 



Tcal 



(5) 



Note here that Vpot is the atomic Coulomb potential, 
while the Vmany-body is the many body potential aver- 
aged from the periodic potential of the atomic arrange- 
ment. The corresponding total energy from Eq. ([5]) is 
given by 



-Eq i C — -Ekin + V'pot + Knany-body 
= Eq ± El + T4nany-body 
^Eo± Ef'K 



(6) 



In this case, Ef'^^ 



is the ionization energy of an atom in 
a many-atom system (not isolated). The exact values of 
Ej are known for an isolated atom. As a consequence, we 
can arrive at Eq. ([2|) from Eq. ([6|) . Apparently, we cannot 
use Eq. ^ to isolate the electronic and phonon contribu- 
tions because we have defined the ^ as a function of the 
Coulomb potential (Vpot), many-body (Knany-body) and 
kinetic (-Bkin) energies. Consequently, the total energy 
can also be rewritten as (from Eq. ([6])) 



cal 



where, j is the sum over the constituent elements in a 
particular compound. For a C nanotube with only one 
type of atom, Eq. ([7]) can be rewritten as 



E^Eo±l3Y,Eu 



(8) 



In Eq. (0), we have defined here that [3 = I + -^^^j 
where (V^(r)) is the averaged many-body potential value. 
Apart from that, the total energy equation for a free- 
electron system is given by 



ircal 

i 



Eq ± [Skin — Ef 



E]iin ^ V^ot ^ ^many — body 
£^kin + Vtotal 

Emu implies free electrons 







pot 



Vr 



many — b 



ody] 



<^ for electrons ± 



(9) 



In Eq. ^ we have substituted Eq. (O for because 
the concept of ionization energy is irrelevant for free- 
electron metals, which do not require excitations from its 
parent atom to conduct electricity. As such, the carrier 
density is independent of temperature and the scattering 
rate is the one that determines the resistivity with respect 
to temperature, impurities, defects, electron-electron and 
electron-plionon interactions. In summary, the total en- 
ergy from Eq. ([2|) carries the fingerprint of each C atom in 
a nanotube and it refers to the difference in the energy 
levels of each atom rather than the absolute values of 
each energy level (eigenvalues) in each atom. Hence, the 
kinetic energy of each electron from each atom will be 
captured by the total energy and preserves the atomic 
level electronic- fingerprint in the nanotube. Using this 
newly defined total energy, we can derive the ionization 
energy based Fermi-Dirac statistics (iFDS) as given be- 
low UM 



feiEo,0 



1 



(,[(Eo+0-E%]/kBT 
^(^0'^) = g[BO_(Bo4)]ABT + l- (10) 

where, Ep is the Fermi level at T = and fc^ is the 
Boltzmann constant. However, substituting the same 
atom in a nanotube gives rise to the influence of many- 
body V{r) and in reality, Ej^"-'' cannot be evaluated from 
Eq. Nevertheless, the i^'J'^"' of an atom or ion in 

a nanotube is proportional to the isolated atom and/or 
ion's Ej as given in Eq. ([S]). It is this property that en- 
ables one to predict the variation of fermionic excitation 
probability in C nanotubes. Therefore, one can employ 



the experimental atomic spectra to estimate, ^ = Ej'^ 



Ej. It is emphasized here that Ej is zero for Boltzmann 
particles. As such, one should not assume that the above 
approximation should give the Boltzmann distribution 
function (BDF) as a classical limit. One can indeed ar- 
rive at BDF by first denying the additional constraint by 
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substituting Ej = 0. Importantly. Eq. ((T0|) is the Fermi- 
Dirac statistics derived specifically for strongly correlated 
matter, where it is not applicable for free-electron system 
(or Fermi gas) as shown in Eq. 



fe{Eo,Ei)N,iE)dE 



Now, before we move on, let us re-examine Eq. ([5]) that 
seems to say nothing about i) the atomic arrangements 
and ii) how to isolate the phonon from electronic con- 
tribution. Firstly, Eq. ([5]) is perfectly applicable for any 
atomic arrangements or crystal structures. The reason is 
that we can incorporate Eqs. ^ and (fTO|) for both non- 
bulk system, namely SNTs as well as for bulk system, 
regardless of its specific crystal structures, since these 
two equations can be normalized by employing the ap- 
propriate density-of-states (DOS). However, for non-bulk 
system of several atoms, including SNTs, we need to in- 
corporate the atomic arrangement explicitly because the 
electronic excitation depends on the number of atoms 
along a certain conducting path (developed here). For 
bulk system with the number of atoms of the order of 
10^^, the effect of different crystal structures do not arise 
because the conducting paths are isotropic and the Ej 
here will and can be dressed accordingly to take this 
structural effect into account [H, [H, [HI, HO]. For 
example, pure diamond and graphite will each have dif- 
ferent valence states and electronic polarizabilities (the 
ability of the valence electrons to excite in a particu- 
lar direction in the presence of electric field), in which, 
these differences are due to the different excitations of 
the valence electrons. These different excitations of the 
valence electrons are the ones that have been captured 
by Eq. (fTO)) through Eq. ([5]). Therefore, in our approach, 
the true DOS and/or atomic arrangements of a partic- 
ular system are unnecessary. The price we pay for this 
is that we cannot calculate the many-body eigenstates 
from Eq. ([2]), but note here that we can indeed prove 
Eq. ^ microscopically for real isolated atoms [31|. In 
other words, our input parameter is the isolated atomic- 
energy-level-difference, or defined here as the ionization 
energy (Ej). The theoretical discussion of how Ej affect 
the polarizability can be found in Ref . [2l| . 



The second issue here is how do we isolate the phonon 
from the electronic counterpart? Basically, we cannot 
and there is no reason to, at least for condensed matter 
that violate free-electron metals and for as long as we do 
not apply this formalism to evaluate thermal conductiv- 
ity. The next question is, how Ej is related to electron- 
phonon interaction in the first place? We will answer 
this shortly. The 1-dimensional (ID) DOS is given by 
Ne{E,m) = [£;-i/2(to*/2)1/2] /^^j^ using e = h^P /2ml 

and k denotes the wave vector. The integral to compute 
carrier density and its solution are given by (after making 
use of Eq. CT) 



- 1/2 


'E%-Ei 


exp 





(11) 



Based on Eq. pTjl . suppose that the system is at tem- 
perature T and it has n number of electrons per unit 
volume. Now, imagine that we reduce the magnitude of 
Ej (small enough that it does not increase n), then the 
only parameter that can change is the effective mass of 
the electron, where m* oc Ej, which in turn implies that 
the electron-phonon coupling (Aohph) has been reduced. 
This same argument with small Ei variations can be ap- 
plied at any reasonable temperatures. However, for free 
electron metals, Aehph is defined as the electron-phonon 
scattering, where electrons and phonons can be treated as 
two different entities that scatter each other. In our ap- 
proach, we do not treat the electrons and phonons, even 
in metallic SNTs as separate entities. In addition, Ej in 
this case has no relation with electron-phonon scatter- 
ing. Switching back to the SNT, the charge (g)-gradient 
along a nanotube's length {I) and the drift velocity [vd) 
of charges can be written as [13] 



dq , dl 

— = n-Kde; — = Vd- 

dl dr 

As such, one can write the current (i) as [32 



dq dq dl 
« = — = — X — = rnrdevd 
dr dl dr 



(12) 



(13) 



Now, the resistance for a single conducting path or 
length, of a SNT is 



7^ 

nd 



V 
i 



1 

i JO 



E dr. 



(14) 



d denotes the tube's diameter and E = electric field. 
We also know that m{dvd/ dr) ~ — eE that eventually 
gives Vd = — eEr/m. Finally, one can arrive at the resis- 
tance of a whole nanotube, as given below 



ird- 



(7rd)2E 



Ah ( 2-Krn, 



= T^d- 

7 

1/2 



7r(iA(13062) exp 



Ei^E%- 



TTdp{El) 

Ei 



(15) 



ksT 

We have substituted, m/ne^Te-e for p{Ej) and the 
electron-electron scattering rate, l/re-e = AT'^. The 
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Tei.ph has been neglected because SNTs are not free- 
electron metals, even the metallic ones. However, for heat 
transport, Tei.ph is not negligible. A is the T-independcnt 
scattering rate constant. The numerical value is obtained 
for T = 300 K. The ID resistivity, p{Ei) for nanotubes 
can be written as 



1/2 

V Ah f 2'Km\ „Q/9 



V kB ) 



exp 



Et 



knT 



(16) 
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are essential. The C in Eq. (|17p represents the Carbon 
atom while z denotes the number of valence electrons 
that can be excited, which will eventually contributes 
to the conductance of CNTs in the presence of applied 
voltage. Apart from N (the number of C atoms along 
a conducting path), the number of valence electron that 
arc excited for conduction in the y direction is not equal 
to the X. Meaning, the strength of the resistance or con- 
ductance in their respective y and x directions of a SNT 
originate from the inequality, Ej{Af) < Ef(Af). The 
subscripts, i = 1, 2, ... z and j adds the C atoms, 1, 
2, and so on continuously along its conducting path or 
length. In the previous work on superconductors and fer- 
roelectrics [11 [H [13, [HI , Eq. ([HD was simply written 
as the average ionization energy of a single ion as given 
in 



Er 



(18) 



The relative magnitude of Ej was then calculated 
based on the percentage of dopant to predict the variation 
of p{T) and dielectric constant. On the contrary, SNTs 
with finite length in nanoscale and with only one type of 
atoms namely, C requires Ej in the form of Eq. (jl7p . 




Carbon atom 



FIG. 1: The intrinsic T-dependence of ID resistivity is y-^/^ 
for T above Ei. For T below Ej, p{Ei) is proportional to 
exp(l/T). There are three different curves for different mag- 
nitudes of Ej. 



The calculated curves from Eq. (jl6p are shown in 
Fig. [T] Interestingly, one of the curve {Ej = 150 K) is 
comparable with the experimental data in Ref. (see 
Figure 2b). The calculation of the total average ioniza- 
tion energy (in the respective y and x directions) can be 
carried out with 



FIG. 2: The arrangement of C atoms in a single-wall Carbon 
nanotube is shown schematically. The resistance is strongly 
influenced by the direction {y, x) in which TZy^x is measured. 
Such observation is due to the inequality, Ej{J\f) < Ef{J\f). 
The defined angles, 6 and 7 can be used to compute the length 
as given in the Eqs. (|19|l and (|20|l . respectively. Uc-c denotes 
the length of covalent bond between two C atoms. 



Figure [5] schematically shows the arrangement of C 
atoms in the y and x directions. Taking Oc-c (0.142 nm) 
as the distance between the two C atoms, one can write 



AG-. 



(17) 



i=2 



1. 



(19) 



Unlike ionic bulk systems, CNTs arc ID systems with 
anisotropic conducting paths, which have covalent bonds. 
Consequently, the following definitions and descriptions 



The Ly and 



arc the lengths along the C— C atom's 
bond. Therefore, the experimentally measurable lengths 
(in real space) can be written as (as a function of L) 
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,y COS{0) = 



cos(0)VaA,-1. 



(20) 



_l + 2cos(7) 



a-c-c 



[l + 2cos(7)]^AA, -1. 
i=2 



(21) 



Recall here that the reason I in Eqs. ([^0]) and (pij) 
are written as functions of L is to take into account the 
higher probability of electrons to conduct along the L. 
The angles, 6 is the chiral angle, while 7 = 90° — 0, which 
are also defined in Fig.[2l The subscript, j = 2 indicates 
the sum starts from the second C atom and so on. The 
chiral vector, Ch is given by [2^ 



Tlx = TTdp{Ei) 
= 7rdA(13062) 



X exp 



3lx 



ac-c[l + 2cos(7)] 



BL 



i]Ei-E° 



7rdA( 13062) exp <j 



(24) 



Equation p4|) accommodates the unit for p{Ej'^), 
which is J7 m~^ (because the unit for ID n from Eq. 
is m"-'^). The length, / varies exponentially as a result 
of Eq. . Figure [3] a) and b) indicate the influence 
of length on resistance via Eq. ((24)) . The • in Fig. [3] a) 
and b) represent the experimental data from de Pablo 
et al. for the nanotubc samples with diameters, d = 1.5 
nm and 1.7 nm respectively. The solid lines are based on 
Eq. 



nai + ma.2, 



where ai and a2 denote the 2D graphcne lattice vec- 
tors, while n and m arc integers. Ch can be related to 

x{e = 0°) and y{e = 30°) with 



Chix) = nai + na2, 
Ch{y) = nai. 

Consequently, Eq. p7p in y and x directions can be 
respectively rewritten as 



Ey[c^+]^Y. 



El 



In 



E^C 



z \ac-cCos{9) 



+ 2cos(7) 



(22) 



(23) 



Now, one can actually substitutes either Eq. (22) or 
Eq. (P51) accordingly into Eq. (fT5)) in order to obtain the 
length-dependent resistance. In addition, we can see that 
both Eqs. (P^ and (P5)) are also determined by the chiral 
vectors. 
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FIG. 3: The length-dependent resistance (72.) based on 
Eq. (|24|l (solid lines) are plotted to evaluate the experimental 
data (•) obtained from Ref. [2^ . 



III. ANALYSIS OF n{l) 

The 7^(Z) of free-electron metals with isotropic distri- 
bution of atoms and electrons can be simply derived as 
TL{1) ~ ^, S denotes the cross section area HI]- How- 
ever, CNTs resistance at 300 K, say in the x direction 
should be written as 



Importantly, the fittings in Fig. [3] a) and b) clearly 
demonstrate that Eq. (|24p gives a reasonable approxima- 
tion. With this model at our disposal, one can utilize the 
fitting parameters namely, 7r(iA(13062) = 37 kfi for d = 
1.5 nm whereas 7rdA(13062) = 300 kn for d = 1.7 nm. 
Therefore, A1.5 = 6.01 x 10^ s-^K'^ and A1.7 = 4.30 x 
10^ s~^K~^. As a result of this, the e-e scattering rate 
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for 1.5 nm and 1.7 nm nanotubes are respectively given 
by re_e = 1.85 x 10"" s and t^_^ = 2.58 x lO^^^ s. 
Eventually, the mean free path, 1^ = vp x T^-e = (8.1 x 
10^)(1.85 X 10-1"*) = 15 nm for d = 1.5 nm, and for d 
= 1.7 nm, /e = 2 nm. Here, the Fermi velocity, vp is ob- 
tained from Ref. [l3|- The other fitting parameter, B for 
d = 1.5 nm and 1.7 nm arc found to be 0.47 and 0.18 re- 
spectively. Throughout this resistance calculations, e-ph 
scattering has been neglected in the usual sense, because 
iFDS have had the electrons dressed with Ej. Meaning, 
the excitation of electrons and holes varies with different 
types of atoms (in this case C), identically with the tra- 
ditional methods discussed by Barnett et al. [Ill, Pere- 
beinos et al. [s^l and Ando [s^. Contrary to iFDS, the 
latter methods utilize the free-electron theory and subse- 
quently the e-ph interaction was determined in order to 
couple it with those free-electrons so as to describe the 
excitation of electrons and holes with different types of 
atoms. Consequently, one can notice that Eq. ([M)) does 
not ignore e-ph interactions in any way. In fact, the exis- 
tence of polaronic effect via Ej has been discussed using 



3JI have 



iFDS [20|. Parallel to this, Perebeinos et al. 
also found strong polaronic effect in SNTs as inevitable. 

The properties of phonons and its influence in CNTs 
specifically and other nanostructur es g enerally have been 
discussed extensively in the Refs. j36l ISTI. [ssI . [sqI . [iol l4ll . 
HH, [4^. Apart from that, Chen et al. j44| pointed out 
the possibility of superconductivity and ferromagnetism 
in SNTs doped by a chain of C atoms. Whereas, Ichida 
et al. 45[ have carried out the necessary analysis on the 
relaxation dynamics of photoexcited states in SNTs us- 
ing femtosecond spectroscopy. They found an interesting 
relationship of which, the e-ph interaction increases with 
decreasing tube diameter. Qualitatively, their result ex- 
plains why for small d (1.5 nm), the B (0.47) determined 
earlier is 2.6x larger than the magnitude of B (0.18), 



which is for large d (1.7 nm). Recall here that B corre- 
sponds to Ej, which is associated to the heavier effective 
mass (polaronic effect). In other words, this polaronic ef- 
fect is due to the interaction between non free-electrons 
and phonons, which enhances the effective mass of the 
charge carriers [20j . On the contrary, for the well known 
e-ph interaction in metals, free-electrons and phonons in- 
teract, that eventually gives rise to e-ph scattering. Hav- 
ing said that, we can now compare our predicted values 
for /e (2 to 15 nm) with the values obtained by consid- 
ering the short optical phonon mean- free-path (10 to 20 
nm, for low bias- voltage and d = 1.5 to 2 nm) that limits 
the electrons mean- free-path Q. 



IV. CONCLUSIONS 

In conclusion, the ionization energy based Fermi- 
Dirac statistics has been employed to derive the length- 
dependent resistance in a single-wall Carbon nanotube. 
It has been shown that such dependence is inevitable 
in a low dimensional and non-free-electron systems at 
nanoscales by using the the recent experimental findings. 
In this paper, it is also highlighted that simple equa- 
tions derived using iFDS are able to capture the transport 
properties of single-wall Carbon nanotubes with reason- 
able accuracy. 
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